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We study dynamical decoupling in a multi-qubit setting, where it is combined with quantum logic 
gates. This is illustrated in terms of computation using Heisenberg interactions only, where global 
decoupling pulses commute with the computation. We derive a rigorous error bound on the trace 
distance or fidelity between the desired computational state and the actual time-evolved state, for 
a system subject to coupling to a bounded-strength bath. The bound is expressed in terms of the 
operator norm of the effective Hamiltonian generating the evolution in the presence of decoupling 
and logic operations. We apply the bound to the case of periodic pulse sequences and find that in 
order maintain a constant trace distance or fidelity, the number of cycles - at fixed pulse interval 
and width - scales in inverse proportion to the square of the number of qubits. This sets a scalability 
limit on the protection of quantum computation using periodic dynamical decoupling. 

PACS numbers: 03.67.-a, 02.70.-c, 03.65.Yz, 89.70.+C 



I. INTRODUCTION 

Quantum information processing harbors enormous 
unleashed potential in the form of efficient algorithms 
for classically intractable tasks Perhaps the largest 
hurdle on the way to a realization of this potential is the 
problem of decohcrcnce, which results when a quantum 
system, such as a quantum computer, interacts with an 
uncontrollable environment [2[ . Decoherence reduces the 
information processing capabilities of quantum comput- 
ers to the point where they can be efficiently simulated 
on a classical computer Q . In spite of dramatic progress 
in the form of a theory of fault tolerant quantum error 
correction (e.g., Q), finding methods for overcoming de- 
coherence that are both efficient and practical remains 
an important challenge. An alternative to quantum er- 
ror correction (QEC) that is substantially less resource- 
intensive is dynamical decoupling (DD) [f|. This method 
does not require feedback or the exponential growth in 
the number of qubits typical of fault tolerant (concate- 
nated) QEC. In DD one applies a succession of short 
pulses to the system, designed to decouple it from the 
environment. This can substantially slow down decoher- 
ence, though not halt it completely, in contrast to the 
promises of fault-tolerant QEC. While initially the gen- 
eral theory of DD was developed under the assumption 
of highly idealized (essentially infinitely fast and strong) 
pulses [y, lDj IS| ) subsequent work relaxed these assump- 
tions, showing that DD can still be beneficial in the pres- 
ence of bounded strength controls [9(. In the simplest 
possible DD protocol, known as "periodic DD" (PDD), 
one applies a certain predetermined sequence over and 
over again. While this protocol typically does not work 



as well as random [lCI I 111 Il2| , recursive-deterministic 
[l3| , or hybrid schemes [l4 ] when finite pulse intervals 
and pulse width are accounted for, it has the advantage 
of simplicity. In this work our purpose is to present a 
rigorous analysis of DD, and in particular to derive error 
bounds on its performance in the periodic (PDD) setting. 

With a few exceptions that belong to the realm of ide- 
alized pulses [HI, Q1I [I?], fTcl or to the paradigm of adia- 
batic quantum computing [191 ]. DD studies have focused 
on preserving quantum information (memory), rather 
than processing it (computation). In order to combine 
computation with DD, Ref. introduced three strate- 
gies: The first strategy requires applying the computa- 
tional operations "stroboscopically" , i.e., at the end of 
each decoupling cycle, where the system is momentarily 
decoherence-free. This is conceptually similar to compu- 
tation over error-correcting codes, where a computational 
gat e is applied at the end of an error-correction cycle 
[201 ]. The disadvantage of this "stroboscopic" approach 
is that, in reality, the computational operations take a 
finite time to implement, so that the system decoheres 
while a computational gate is being applied to it. The 
second strategy, is to alternate and modulate the control 
Hamiltonian used to implement quantum computation, 
in which the net overall effect of the DD operations still 
allows a desired unitary operation on the system, along 
with the correction of errors. The third strategy pro- 
posed in Ref. [u| is to use DD pulses that commute 
with the computational operations, so that the two can 
be executed simultaneously. Here we address the prob- 
lem of circuit-model quantum computation [l| using DD 
with realistic pulse assumptions. We combine DD with 
computation via the use of codes and universality results 
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arising from the theory of decoherence-free subspaces and 
subsystems (DFSs) [H, [H, [H . Our method is concep- 
tually related to a hybrid version of the second and third 
strategies of Ref. [15[ , in that we impose the commutation 
condition between the DD pulses and the computational 
Hamiltonian, but we find that improved performance is 
obtained if the DD and computational operations simply 
alternate. Thus, in our scheme the computational gate 
is "spread" over an entire DD cycle (or, conversely, the 
DD cycle is spread over the computational gate). We 
fully incorporate finite pulse intervals and pulse widths 
and assess the performance of our scheme in the PDD 
setting. We find a rigorous error bound, from which it 
follows that for a fixed error the number of DD cycles 
cannot scale faster than the inverse square of the system 
size (at fixed pulse width and pulse interval) . This means 
that there is a tradeoff between the length of time over 
which decoherence errors can be suppressed using PDD, 
and the scalability of a quantum computation it is meant 
to protect. 

The structure of this paper is as follows. We define the 
model in Section [ill We provide background on dynam- 
ical decoupling in Section IIII1 where we also derive the 
effective Hamiltonian describing the evolution under the 
action of decoupling and computation. This leads to the 
"error phase" , namely the effective Hamiltonian times 
time (a type of action), which is the quantity we wish to 
minimize. In Section HVl we derive rigorous error bounds 
that relate the error between the desired and actual final 
state to the norm of the error phase. In Section [V] we 
estimate the error associated with the decoupled evolu- 
tion (i.e., the evolution in the presence of a DD pulse 
sequence), relative to the decoherence-free evolution (no 
system-bath coupling) . Section IVII is where we derive 
our key result: we apply the idea of encoded operations 
and dynamical decoupling to PDD, and compute the er- 
ror bound. In Section [VIII we illustrate our construction 
with encoded DD-computation in a quantum dots set- 
ting, where computation is implemented via Heisenberg 
interactions. We conclude with a discussion of our re- 
sults in Section IVIIII Extensive background material is 
presented in the appendices. 



II. MODEL 

We express the total Hamiltonian for system plus bath 
in the form 



H(t) = H ctll {t) (g>I B + H CII + I S ®H B 



(1) 



where / is the identity operator, H ctv i acts on the system 
only and serves to implement (encoded) control opera- 
tions such as logic gates, H crr is the "error" Hamiltonian 
(system-bath couplings, undesired interactions among 
system qubits that do not commute with H ctT i), and Hb 
is the pure-bath Hamiltonian. Let f/ c tri be the (encoded) 
logic gate generated by switching on H ctT \ for duration T, 



in the absence of the bath and any undesired interactions 
within the system: 



U ctI i{T) = Texp 



-1 I H ctl \{t)dt 




= exp(-i9R), (2) 



where T denotes time-ordering, R is a dimensionless logic 
operator, and 9 is the angle of rotation around this op- 
erator. However, due to the presence of the undesired 
-ffcrr and Hb terms, we will in fact obtain the follow- 
ing unitary acting on the joint system and bath Hilbcrt 
space: 



[^barcCr) = Texp 



(H cM (t) + H e „ + H B ) dt 



This is the essence of the problem of any quantum con- 
trol procedure, whether it be for quantum information 
processing or other purposes: [/bare entangles system and 
bath and implements a tranformation on the system that 
can be very different from the desired U c t T \. Our goal in 
this work is to show how to modify [/bare so that the 
distance between a state evolving under it and a state 
evolving under f7 c tri can be made arbitrarily small. This 
will be done by adding another Hamiltonian to the sys- 
tem, which implement DD operations, and is designed to 
effectively cancel H crT without interfering with H ctl i. 



III. DYNAMICAL DECOUPLING 
BACKGROUND AND THE ERROR PHASE 

A. Dynamical Decoupling Denned 

We assume that the decoupling operations are realized 
as pulses Pi by a time-dependent Hamiltonian Hw(t). 
The essential condition that will ensure that the decou- 
pling pulses interfere minimally with the control opera- 
tions is: 



[#DD(*),flctrl(t')]=0 W,t'. 



(4) 



The total propagator is now generated by the time- 
dependent total Hamiltonian 



H tot (t) = F DD (t) + #ctri (t) + H eII + H B , (5) 



i.e., 



U(T) =Texp 



H tot (t)dt 



The pulses are applied at times {tj}f_ given by 



j=0 



U 



j(r + S), 



(6) 



(7) 



where r is the pulse interval and 5 is the pulse width. 
From hereon we assume for simplicity that Huuit) is 
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piecewise constant (thus we are performing a worst-case 
analysis: our conclusions can only be improved by pulse 
shaping [H [H): 



.? e {0,---,iV-i}. 



Let 



T = % = N(t + 5) 



(9) 



denote the time it takes to complete one DD cycle, con- 
sisting of N pulses generated by the H 3 p : 



P 3 = exp(-iH^5) j = 1, N. 



(10) 



The commutation condition (|4|) becomes 

[H^,H ctll (t)\ = [Pj,H ctTl (t)} = Vj. (11) 

This will allow us to import many of the results of the 
control-free scenario, i.e., when H ctI \(t) = 0. For the 
remainder of this section we review this setting, and re- 
turn to the question of how to ensure the commutation 
condition in section IVTTl 

Denoting a free evolution period (when i^DD = 0) of 
duration r by 



f T = exp [-ir (H eTI + H B )} , 
a single cycle can be written as 



(12) 



P N f T P N -l£ tPn-2 ■ ■ ■ Plfr = P N fr[PlPN}PN-lfr[{PNPN-l)HPNPN-l)}PN-2 ■ ■ • P%f r 

= (DrrfrDlfKDff^frD^jDjr-i ■ ■ ■ £>if T I>1 



(13) 



where the unitary "decoupling group" @ Q = {D j}jLi 
has elements defined as 



II /'v •••/'.,• D 1 =I S , (14) 

where the condition D\ = Ig is imposed because of the 
appearance of D[ in Eq. (| 13|) : this imposes a relation 
among the pulse Hamiltonians Hp via Eq. (jTU|) . Note 
that this is only possible in the zero width limit, since 
such a relation cannot be satisfied when the system-bath 
and bath Hamiltonians are present during the pulse. 

The effective Hamiltonian ffff (T) can be computed 
using the Magnus expansion [26[ (see also Appendix [A"|l . 
To first order in the Magnus expansion 

N 

HesiT) = ^D^H^ + H^D] 

= Il g (H eII ) + NH B . (15) 
This can be viewed as a projection 

N 

Ug{H exT ) = Y,D j H exT D] (16) 
into the centralizer Z(Q) = {V\ [V,Dj] = 0, VDj E Q}. 



Indeed, since [H B ,Dj] — for all j, 

N 

DtWi$,D k ] = D\ ^2[DjH eII Dt + H B , D k ] 

N 

= Y,( D h D A H °"( D h D tf - D i H *" D ) 

i=i 
= 

=> [H^,D k ]=0 Vfc, (17) 

where in the last equality we used the group closure prop- 
erty: Vfc, j 3i such that DuDj = Di. For a unitary ir- 
reducible representation of Q this immediately implies, 
by Schur's lemma, that HQ = ds ® H B where c is a 
constant, i.e., H\) acts harmlessly on the system. 

If the algebraic structure of H e „ is known, we can 
choose the decoupling group such that it satisfies the "de- 
coupling condition" for any Hamiltonian H elT : [1,0] 

N 

He; (fferr) = ^ Pjll.-D, = 0. (18) 

i=i 

In the limit r, S — > 0, and in the absence of control, the 
first order Magnus expansion is exact and condition [)18p 
guarantees the stroboscopic elimination of H m , in the 
sense that H^(T) — 0, and this would be true at the 
end of every DD cycle. Another way to understand con- 
dition (fT5|) is to recall that D\ = Is, which means that 
J2j=2 DjHenD'j = —H elt : the negative sign in front of 
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flerr means that the role of the decoupling group is to 
effectively time-reverse the error Hamiltonian at the end 
of the cycle. 



B. Interaction Picture 

In a setting where decoupling works perfectly the sys- 
tem evolves independently from the bath, purely under 
the action of the control Hamiltonian. Therefore we use 
the interaction picture of 



H scc = -f/ctri + Hb 



(19) 



(the sum of the secular terms) to calculate the full prop- 
agator [Eq. ((6|)]: 



U(t) = u scc u CTT (t,o), 



where 



sec — 

t 



^ctri(i) = Texp 



— i 



H ctv \(s)ds 



U B {t) = exp(-itH B ). 



(20) 

(21) 
(22) 
(23) 



If |*(<)) = U(t)\V(0)) (Schrodinger picture) then 
= U^ cc \V(t)) = U eII (t)\V(0)) is the correspond- 
ing state in the interaction picture. Similarly, for mixed 
states: p(t) = U^ cc p(t)U scc = U eTT (t, 0)p(0) 0). The 
interaction picture propagator C/ crr contains all the "er- 
rors", in the sense that if it becomes the identity op- 
erator then decoupling is perfect. For then U{t) — 
U c tii{t) <8> Unit) and the desired system dynamics is com- 
pletely decoupled from the bath. In this sense the in- 
teraction picture is naturally suited to our analysis: by 
moving the "ideal" evolution U c t r \(t) ®> Unit) to the left, 
we have isolated the "error propagator" U CII . 
C/ crr satisfies the Schrodinger equation 



dU e „(t, 0) 



dt 



= -%H m (t)U m (t,0), U eiT (0,0) =1, (24) 



with 



H m (t) = U B Ul M [H m (t) + H eiI ]U B U cttl 
= H BB (t) + Ad tHs jH m ], 



(25) 



where the linear adjoint map Ad^[i3] has the Baker- 
Campbell-Hausdorff formula |27| 



Ad A [B] ee e~ lA Be %A = 



(-»)* 



(26) 



where [ n A, B] denotes a nested commutator term 
[A, [■ ■ ■ [A, B] in which A appears n times. Note that 
- thanks to the commutation condition - Huoit) re- 
mains invariant under the interaction picture transfor- 
mation in Eq. (|25|) . This is where the commutation con- 
dition shows up explicitly in our analysis. 



Let us define an effective (dimensionless) "error phase" 
via [281 : 



expH$ E (T)] ee U e „(T,0). 



(27) 



Thus $e(T) is the final effective Hamiltonian times the 
total time, and it measures the deviation from ideal dy- 
namics. In other words, the goal of the decoupling pro- 
cedure is to minimize $b(T). In Section HVl we relate $e 
to conventional fidelity measures. Throughout this work 
we repeatedly use the technique of expressing unitaries 
in terms of the "final effective Hamiltonian" . In fact, this 
was already done in our review of DD above, when we 
used the effective Hamiltonian H^g (T) in Eq. (| 13|) . 



C. The Error Phase 

We now wish to calculate the total propagator U 
[Eq. ©] in the presence of both decoupling and con- 
trol. The evolution generated by H eir (t) [Eq. ()25f)] can 
be decomposed into "free" and pulse periods as follows: 



H e „(t) 

(Ad tHsoc (ffcrr) U-X < t < ti-x + T 

\ H® + Ad tHacc (H m ) t 4 _i + t < t < U 



(28) 



In Appendix [B] we prove the following "switching 
lemma" : 

Lemma 1 The propagator generated by a "switched 
Hamiltonian" 

H(t) = Hi(t) U-i<t<U, i = l,...,N, (29) 

can be decomposed into corresponding segments: 

U(t N ,t ) = U(t N ,t If -i)--'U(t 1 ,to), (30) 

where U(ti+i,ti), with U < t < tj+i, satisfies the 
Schrodinger equation 



dU(t, U) 
dt 



= -iH{t)U{t,U), U{U,U) = I. (31) 



We can thus write: 

U eTI {T, 0) = U CII {t N ,t N ^) ■ ■ • [/ err (ti,t ), (32) 

where U eII (ti, ti-i), with tj_i < t < ti, satisfies the 
Schrodinger equation 



dU crr (t,tj-i) 
dt 



= -iff err (t)f/err(Mi-l) (33) 



Uen{ti-1, U-i) — I, 

whose formal solution is 



U m (t,t i - 1 )=Texp[-i H m (t)dt]. (34) 

Jti- X 
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Analogously, we can further decompose each segment 
into a pulse and a free evolution, each with an effective 
Hamiltonian. The pulse part is 

U eTT (ti, ti-i) = U eTT (ti,ti — S)U eTI (ti — 5,ti-i) 

= Te'*^-* *" At)dt U m (U - S, ti-i) 

= P i e- WH «I/ Brr (t i -*,t i _i), (35) 

which, using Eq. (|28|) . serves to define the effective error 

expl-iSHgiS)] ee P }Te- i %-^ )+Ad ^ H °"» dt 

(36) 

associated with the width of the ideal pulse P{ — 
exp(-i8H^). 

The i-th free segment (£4—1, U~5) is similarly generated 

(i) 

by an effective Hamiltonian HcA defined via 

exp[-iTH^l] = ?7crr(£i-l +T,U-l) 

— 7" e _i /o T Ad (»+«i-l)H S ac[ ff «']^ s _ (37) 

The overall error unitary U CII (T,0) can thus be written 
as 

U eII (T,0) = P N ew[-i5Hg]exp[-iTH<g)] 

■•• x Piexp[-iSH^]exp[-iTH^}] (38) 

To incorporate the effect of the DD operations we recall 
the definition of the decoupling group in terms of the 
pulse unitaries [Eq. |[T3|)] and rewrite Eq. (f3"8")) as 

C/ crr (T,0) = D N e- i5H °" DlfDjfe-^^ Dlf 
... x D^-^DlD^-^Dl 

N 

= \{ C M-tiD]H^D 3 ). (39) 

3=1 

By Lemma [T], the following time-dependent Hamiltonian 
generates U e „(T, 0): 

[DJ&D\ iorU<t<U + r 
[DiH^Dl ioYU+T<t< t l+1 

The "free evolution crror-Hamiltonian" H^A and "pulse 
error-Hamiltonian" arc defined, respectively, in Eqs. 
(I5T1) and (J3SJ) • Gathering our results we can write: 

cxp H$ £ (T)] = U m (T,0) 

= TcxpH / H m (t)dt}. (41) 
Jo 

IV. ERROR BOUNDS 

In this section we derive rigorous error bounds that 
relate the error between the desired and actual final state 



to the norm of the error phase $>e(T). Throughout this 
work we use the trace distance 

D( Pl ,p 2 ) = ^\\ Pl - p 2 \\ l7 (42) 

where 

||i4||i=tr(VAtA), (43) 

as the distance measure between state, and the quantum 
fidelity, defined for any pair of positive operators A and 
B: 

Fq (A, B) = WyfAsfBh AW = st=s tiV VBAVB. (44) 

There is a useful relation between the trace distance and 
the quantum fidelity [l[ : 

1 -D(p u p 2 ) <F Q {p u p 2 )< Vl-£>(pi,/>2) 2 , (45) 

which means that the trace distance and fidelity can be 
used to bound one another from below and above. 

When one or more of the states is pure (|^>i), ^2)), 
we shall write D(ipi,ip2) and F(ipi,ip2), or use a mixed 
notation D(tpi,p^) and F(ipi, P 2), etc. 

We also make repeated use of the operator norm 

PIU- sup \\A\1>)\\. (46) 

11-011 = 1 

For a review of these measures along with key properties 
see Appendix ICl 

In the absence of the bath the control Hamiltonian 
H c t T \(t) would implement a quantum computation via 
the propagator E/ C tri [Eq. Equivalently, the state 

of the quantum computer at the final time T would 
be described by the solution \ip(T)) of the Schrodinger 
equation = — iH ctl i\ip). Imperfect control of if ctr i(t) 
means that even in the absence of the bath, \ip(T)) is not 
the ideal final state, which would be obtained if one could 
implement a completely accurate and noise-free Hamil- 
tonian i?ctri al (0; with corresponding final state \<p(T)). 
Minimization of the corresponding closed-system control 
error 

S id = D[^(T),0(T)] (47) 

belongs to the realm of fault-tolerant quantum compu- 
tation Q and composite pulse techniques [3(|, and will 
not be addressed here. 

The initial bath state is pB (0) and in the absence of 
coupling to the system it evolves under the pure-bath 
Hamiltonian to p° B (t;9) = U B {t)p B {9)U B {t)^ ( tne super- 
script denotes no system-bath coupling). 

In the general mixed-state setting we distinguish be- 
tween "ideal" system evolution described by a pure state 

P ^(£) = m))w% .^th \m = u^mm and 

^ctri al generated by ^tri ' an< ^ bath- free non-ideal sys- 
tem evolution (due to control errors), described by a 
mixed state p%(t) (the mixed nature can be due to, e.g., 



6 



the need to average over stochastic realizations of uni- 
tary evolutions). In the absence of any coupling between 
system and bath the joint initial states Pg eal (0) ® Ps(0) 
or ^5(0)0^^(0) evolve in the two scenarios to /9 ldeal (i) = 
pf eal (t) ® p° B (t) or p°(t) = p° s {t) <gi respectively. 

Then the final error due to imperfect control in the 
uncoupled setting is 

Ad = D[p°(T),p idca \T)] = D[p|(T),^ deal (T)], (48) 



where we have used the multiplicativity property (|C7j) 
and unitary invariance. Minimization of the pure-system 
control error Ad [generalization of Eq. (|47[) ] once again 
lies in the domain of fault tolerant quantum error correc- 
tion [5^1 and composite pulse techniques [30| . 

The actual system-bath state obtained by time evolu- 
tion under the full propagator U(t) [Eq. ©] is p(t) — 
U(t)[ps(0) <8> pB(0)]U(ty , and the actual system state is 
ps(t) = tr Bp(t). The distance we wish to minimize is the 
distance between the actual final system state p{t) and the 
ideal final system state (no control errors, no coupling to 
the bath): 



A 



D[p s (n P r al (T) 



Define 



D tot = D[p{T) lP ^\T)} 



(49) 



(50) 



By virtue of Eq. (jClip we know that removing the partial 
trace can only increase the distance between states, i.e., 



As < A 



(51) 



Let 



Add ee D[p(T), p°(T)} = D[p(T), p»], (52) 



where we have used the fact that in the interaction pic- 
ture p°(t) = pP(0) = p s (0) ® pb(0). Add is the distance 
due to coupling between system and bath, and the role 
of the decoupling procedure is to minimize this distance. 
Using the triangle inequality on \\p(T) — p°(T)+p°(T) — 



idcal 



(T)||i we have 



Aot < Ddd + Ad, 



(53) 



which shows that minimizing the total error can be done 
by separately minimizing the open-system decoupling er- 
ror and the closed-system control error. 

In Appendix [D] (see also [3l| for a more general treat- 
ment) we prove the following lemma: 

Lemma 2 Let U = exp(— iA) where A is hermitian. 
Then for any submultiplicative norm 

\\UBU ] -B\\ < \\B\\ min[2,e 2||A|l °° - 1] (54) 

2||-4||oo <1 

< 2||B||min[l,(e-l)||A|| 00 ] 



By identifying A with the error phase $>e{T) and B 
with po(0) this allows us to write 



D 



DD 



i||f/ OI ,.(r,0)p(0)C/t rr (T,0) - po(0)||i 
< illpoCOJHimmp, ( e 2|l*-C^)H- _ 1)] 



min[l, I( e 2||i» S (r)|U _!)] 



2||* B (T)|U<1 

< 2||$ B (T)|| 



(55) 



Inequality (|55[) shows that minimization of the error 
phase <& e (T) is sufficient for minimizing the decoherence 
error Add- Combining our bounds [Eqs. (gSJ), (|55|) . and 
(|55p] we have the quantum fidelity lower bound between 
the actual and ideal system state: 



F Q [p s (T), P r ai (T)} 



(56) 



> 1 - DlpUnp't^T)} min[l, -(e 2 ^*™~ 1)]. 

In other words, minimization of the pure-system control 
distance together with minimization of the error phase 
$ e (T) is sufficient for minimization of the total distance 
As- Note that the bound we have derived is not neces- 
sarily tight: it is possible to minimize Add and Ad si- 
multaneously, rather than separately, as is done in fault 
tolerant quantum error correction [29j |. 



V. ERROR ESTIMATES FOR DYNAMICALLY 
DECOUPLED LOGIC GATES 

Our goal in this section is to estimate the error as- 
sociated with the decoupled evolution (i.e., the evolu- 
tion in the presence of a DD pulse sequence), relative 
to the decoherence-free evolution (no system-bath cou- 
pling) . The decoupled evolution at the end of a DD cycle 
is described by the propagator U crr (T, 0) of Eq. (f3"9"| . The 
appropriate dimensionless error parameter is the norm of 
the total error phase ||$e||oo [Eq- ([27)) ]. Our strategy for 
estimating $ e will be to calculate approximations to the 
final effective Hamiltonian and then to bound its norm. 
As we showed in Section ITVl in the limit that $ e vanishes 
the final state is free of decoherence errors. 

Our main technical tool is the following lemma. For a 
proof see [31 1. 



Lemma 3 Consider a quantum evolution generated by 
a time- dependent Hamiltonian H(s) = Hq(s) + V(s), 
< s < t, with propagators satisfying dU(s,0)/ds = 
-iH(s)U(s,0) and dU (s,0) = -iH (s)U (s,0). Then 
there exists a time- dependent Hamiltonian H e ff(t) such 
that 



eM-HH eff (t)] ee U^{t,0)U(t,0). 



(57) 
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and the following inequality holds for any unitarily in- 
variant norm: 



\Heff\\ 



1 

< - 



ds\\V(s)\\ = (\\V\\) t 



< sup ||V(s)||. 

0<s<t 



(58) 
(59) 



This lemma allows us to relate the strength of the ef- 
fective interaction picture Hamiltonian at the end of the 
evolution H g(t) to the strength of the (time-dependent) 
perturbation V. 

As a hrst application, let us relate ||-H^.||<x3 to 
1 1 -Herri | cx>- Comparing lemma [3] with Eq. (|36| and iden- 
tifying Ho with Hp [and thus Uo{t, 0) with the ideal 
pulse Pi], V(t) with AdtH BBC (H eiT ) [and thus H(t) with 
H^+AdtH^iH^)], and 5H eS (6) with 5H^ T (S), we have 

\\H*(S)\\«> < (\\Horr\Us < 



sup \\H erl (t)\\ 

t i+1 -8<t<t i+1 



, . (60) 

This means that the application of a pulse, with inclu- 
sion of the system-bath coupling during the pulse as in 
Eq. (l36|) . does not cause a growth in the error rate. This 
is rather remarkable and can be summarized as "pulses 
can't hurt". 1 

Similarly, by setting H = and V(t) =AdtH eec (H elI ) 
in lemma [31 we obtain for the free evolution: 

Halloo < <ll#err||oc>T < Sup || H e „(t) || (61) 

ti<t<ti+T 

Now let us return to H m [Eq. (|40f ], i.e., the Hamil- 
tonian describing the total evolution over a DD cy- 
cle. From now on we simply denote (||-H e rr||ao)i5 an d 
( 1 1 -H crr 1 1 oo ) T by ||i?crr||oo- Then the last two inequali- 



ties yield HAJT&^illoo < \\H e „\\ x and || D.HS'DjWoo < 



W nti 



-Herri loo, SO that 



\H m (t)\\oo < ll^err) 



(62) 



At this point we are ready to use the Magnus expansion 
to estimate the error phase $^(T). Recalling Eq. (|4"Tj) . 
the Magnus expansion for the error phase is given by 

®e(T) = f d Sl H m { Sl ) (63) 
Jo 

+ o / ds i [ ds 2 [H m (s 2 ),H m (s 1 )]-\ , 

z Jo Jo 

and converges as long as J„ T dsi||-H m (si)|| 00 < 7r 
[Eq. (|A7|) ]. i.e., a sufheient condition for convergence is 

ryifetrlloo < 7T. (64) 



1 Of course by other measures pulses can hurt. For example, 
broadband pulses can cause unwanted transitions. But for the 
purposes of our error estimates the absence of growth of the error 
norm is the crucial aspect. 



We assume that our decoupling sequence {Pi} is de- 
signed for cancelling error terms up to the first order in 
the Magnus expansion, as in Section HTT1 Accordingly we 
rewrite &e as a sum of the first order terms and a second 
order correction, which we can in principle improve upon 
by designing a pulse sequence that cancels error terms 
up to a higher order: 



*b(T) 



f H m { S )ds + $2nd(T). (65) 

Jo 



In Appendix [E] we prove the following lemma 

Lemma 4 Consider a time- dependent Hamiltonian 
H(t), < t < T, and the partial sum of kth and higher 
order terms in the corresponding Magnus expansion: 



$1 



i—k 



Assume the Magnus expansion converges. Then 
Halloo < c k (T sup ||!Z(t)||oo) fc 

0<t<T 

where Ck = 0(1) is a constant. 
Thus, subject to Eq. 

||*2nd(T)||oo < c(T||fletr||oo) 2 



(66) 



(67) 



(68) 



for some constant c. The fact that starting from ar- 
bitrary fcth order the error phase is upper bounded 
by (T||i/ crr || oc ) fe means that the design of higher or- 
der pulse sequences can be very advantageous in achiev- 
ing improved convergence of the DD procedure (see also 
[25I I32T l37l|). but we will not pursue this here. 

To calculate the first order integral in Eq. (|65f we sep- 
arate the pulse and free parts: 



/ H m (s)ds = $pulso + $free, 

Jo 

/ H m (s)ds, 
'u-s 

rU- 

E 



pulse 



N 

N-l rti+T 



free 



-n Jt 



i=0 



H m (s)ds. 



(69) 



Using Eq . (f6"2"]) . the error phase due to the pulses ^puise is 
bounded by 



where 



^pulsclloo — ^||-H C 



NS 



(70) 



(71) 



is the total length of the pulse durations. Without use 
of additional techniques such as pulse shaping [25|, [HJ 
or composite pulse sequences [30j , the only means at our 
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disposal to minimize this error is to make the pulse du- 
ration S small. Taking stock, we have, so far: 

HSuCniloo < || I ffmMdslloo + II^COHoo (72) 

Jo 

< AUfferrlloo + ||$f r co||oo + c(T||F er r||oo) 2 . 

The "free error" H^froelloo is the target of the DD 
pulses. Explicitly, we have, using Eq. (j^Djl : 



JV 



(73) 



i=l 



The effective error Hamiltonians Hell can be Magnus- 
expanded to first order in H crl (t) =AdtH Bec [H e „] [recall 
Eq. ([28]) ]. so that the higher order commutators arising 
from the time-ordering in the definition of Hell, are in- 
cluded in a new term C that will be absorbed into $2nd- 
First, 

exp[-irH^\ = r expH / H e „(t)dt] 

rti-i+T 

= expH / H c „(t)dt + C w ] 

r t«-l+T oo 



cxp [ — I 



E 



1 n=0 



[ Ti iT sec ,# err ] < fi + C'«], (74) 

where by Lemma |4] 

Halloo < d(T Sup ||i/ crr (t)||oo) 2 

tj-l<t<ti-l+T 

- rfCrH^errlloc) 2 , (75) 

and where d is a constant. The integrals yield 
/^- 1+T (-it)" 

Jtj-x 

(-i) T 



-rff 



r(i) 



(76) 



Therefore we can define ifcrr via the expansion 

oo 

r/rW = Tif err + 2 [»# 8 ec , #err] + C W • (77) 
n=l 

Returning now to $f lcc [Eq. (T73"j) ], notice that the first 
term (H erl ) in Eq.([77j) docs not depend on fj and is sin- 
gled out, so that we may write: 

$free = ^dec + $undcc + C, 
N 



$ dcc = T^DiH eiI Dl 

i=l 
N oo 

^undcc — ^ ^ fn.i[nH SC c: ^iH e rr^i] 



i—1 n—1 

N 



c = ^C w , 



(78) 



where 



||C|L < Nd(T\\H cl . r \U 2 
1 



N 



[d(T- A)||fl wr || 0Q )'' 



(79) 



The purpose of the DD procedure is, of course, to cancel 
<&dcc [recall Eq. (TTS)) ] . Pulse sequences that cancel higher 
order terms (n > 1 in Dundee) can be found, but this will 
not be pursued here. Thus in our case the undecoupled 
terms will be given by $ un d C c- Define 



P=\\H S , 



J=\\He 



(80) 



Let us first note that, using the triangle inequality on 
Dundee = $free - $dcc and Eqs. $73J) , f77|) : 

PuBdeclU < PfreelL 
N 

< r^HA^^llloo 

N 



r^||i/W||oo<T||i/ crl .|| 



= JT. 



(81) 



This trivial upper bound simply means that the undecou- 
pled error is bounded above by "do nothing". Another 
upper bound for Dundee can be found by making use of 
norm submultiplicativity: 



U,B)\\ 00 <2\\A\\ 00 \\B^ 



(82) 



for any pair of operators A and B in the combined 
system-bath Hilbert space. Then: 



|$undec||oo < EE'^' 1 1 [n^scc , A^crr ^ 1 1 



i—1 n—1 

N oo 



< 52J2\fn,i\m n J 



i—1 n—1 

N c 



i—1 n—1 

N oo 



(n + l)\ 



^ J EE ^xrv — — ( w 



j ( 2 /^) xn+l _ y.n+1 

J Z 1(n + 1)! 2_ r t i 

n—1 v 7 i=l 

</ ^n+1 

2^E (^+T)T( 2 ^ n+1 

n—1 v 7 



AT 



JT 



2/3T 



(83) 



where in the penultimate equality we used T = tjy and 
< = 0. 
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Combining the bounds in Eqs. (|8Tj) . (|83|) we obtain 
the following bound on the strength of the undecoupled 
terms: 



undo 



< JT min 



exp(2ffT)-l 
2/3T 



(84) 



Combining the expressions for various parts of the total 
error phase $£,we obtain the following upper bound: 



1$; 



< ||$2nd||oc + Impulse Hoc 4 

dec||oo H - ||C||oo 

< c(JT) 2 + JA + 
exp(2/3T) 



+JT rain 



1 



1, 



2/3T 



^[dJ(T-A)) 2 ]. 



(85) 



Now note that for fixed N and A we can always write 
d(T — A)/\/N as c'T, where c' accounts for the shift and 
rescaling of T. This allows us to absorb jf[d(T — A) J) 2 ] 
into c{TJ) 2 (redefining c in the process), so that: 



exp(2ffT)-l 
2/3T 



|«b|L < c(JT) 2 + JA+JT min 

z ■ > ; 

This, in conjunction with Eq. (|56[) . finally gives us the 
desired lower bound on the quantum fidelity of one period 
of DD. 



VI. PERIODIC DYNAMICAL DECOUPLING 



Note that in principle the bound (|86p is appropri- 
ate for any DD sequence, since the time T is arbitrary 
(subject to the convergence of the Magnus expansion) 
and the decoupling group can have arbitrarily many el- 
ements. However, in practice DD pulse sequences have 
some deterministic structure, such as periodicity or self- 
similarity, or are random. Structure generally results in 
im p rov ed p erformance under appropriate circumstances 
[T3, [H, [3|3E HI, H3, HI] , and hence the bound (JSBJ) may 
be too weak. 

In this section we apply the idea of encoded operations 
and dynamical decoupling to the periodic case (PDD) 
and derive the final-time error bound. The encoded oper- 
ation consists of the switching of a physical Hamiltonian 
corresponding to a logical Hamiltonian for a duration of 
T m . This switching period T m is punctuated at various 
points by the action of dynamical decoupling operations. 
In the previous section, the analysis was performed for a 
basic cycle of N pulses. In this section we consider what 
happens when this sequence is applied m times. 

Consider a basic decoupling sequence p designed to 
cancel all terms in H erl , as in Eq. (fT5)) : 



p (f T ) = PNfrPN-xfrPN-2 • ■ • Pit r, 



(87) 



where {Pi} is the sequence of N decoupling pulses and 
f T denotes a "pause" of duration r = in decoupling, 
during which the control Hamiltonian H ctr \(t) generat- 
ing the encoded logic gate is operative. Consider now 
the longer periodic sequence PDD m formed by repeating 
p (f r ) m times to obtain a sequence of length Ti ong = mT 
with N m — mN pulses: 



PDD m = n^( f -) 



(88) 



In the absence of encoded operations the sequence p (f T ) 
is designed to cancel dynamics up to the first order. The 
longer sequence PDD m has the same canceling properties 
as the sequence p in the limit of r — * 0. 

So far we have not been specific about how we imple- 
ment the encoded operation. Namely, we have consid- 
ered general time-dependent control Hamiltonians. For 
simplicity, from now on we consider the following simple 
method for realizing encoded operations. First, we only 
implement one logic gate during each PDD sequence. 
In other words, a new logic gate requires a new PDD 
sequence. Second, each logic gate is implemented in 
terms of a constant control Hamiltonian. Thus, if ideally 
we wish to implement U c t r i(T m ) = exp(— iT m _ff ctr i) = 
exp(—i8R) [Eq. @], where H ctl \ = XR with A the mag- 
nitude of i/ctri and 9 = XT m the phase, then in practice 
we will implement the decoupling-free intervals as 



f T = exp 



N m 



(89) 



I.e., the encoded operation is implemented little- by- little, 
using 7Vi on g equal iV m th root segments. 

Let us now find a bound on the fidelity of PDD m in 
this setting. Since we implement the encoded operations 
using the fixed step f T , the propagator for each cycle 
in the periodic sequence is the same, and hence so is 
the error phase at the end of each DD cycle. Formally, 
the total propagator in the interaction picture is simply 
[recall Eq. (|4T|)]: 



C/ O rr(T m ,0) = U m (T m ,(m-l)T) 

■■■U eiI {2T,T)U e „{T,0) 



i<S> E UT) 



3=1 
-im*a(T) 



-i* E (T)lra 



-i<E> F 



(90) 



Recalling our fidelity bound Eq. (|56j) . our task is to 
estimate the norm of the error phase associated with the 
periodic sequence after time T m , i.e., $pdd,„ = w$^;(T). 
It thus follows immediately from Eq. (|86|) that: 



|$pddJ|oo < c(JT m ) 2 /m + N m JS 
+JT m min 



(91) 



1 exp(2/3T m /m) - 1 _ . 



2p(T m /m) 
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In the limit of S = and f3T m < 1, we have (second 
order Taylor expansion): 

||$pddJ|oc < m(cJ 2 + J0)T 2 . (92) 

We postpone an analysis of this result until Sec- 
tion Enn 

VII. EXAMPLE: QUANTUM COMPUTATION 
USING THE HEISENBERG INTERACTION 

The commutation condition Q is crucial to our re- 
sults. At first sight it appears that one cannot satisfy it 
while having non-trivial decoupling operations. However, 
as pointed out in Ref. |l9j, it can be satisfied using the 
double commutant construction, which we now explain. 

The decoupling group Q induces a decomposition of 
the system Hilbert space Hs via i ts g roup algebra CQ 
and its commutant CQ' , as follows [33. l40||: 

H s = 0C" J ®C dj , (93) 
./ 

C $ - ® M d j > Cg' S M„ J <X> J dj . (94) 
j j 

Here nj and are, respectively, the multiplicity and 
dimension of the Jth irrep of the unitary representation 
chosen for Q, while In and Mjv are, respectively, the N x 
N identity matrix and unspecified complex- valued N x N 
matrices. We encode the computational state into (one 
of) the left factors Cj = C nj , i.e., each such factor (with 
J fixed) represents an nj-dimensional code Cj storing 
log d nj quofits. Our DD pulses act on the right factors. 
As shown in [3§| , the dynamically decoupled evolution on 
each factor (code) Cj will be noiseless in the ideal limit 
w,t -> iff Ug(S a ) = 7 Xj, a I nj ® Idj [the projection 
was defined in Eq. (fT()|) ] for all system operators S a in 

H S b, whence H£ = 7 [(I nj ® idj)]^® E Q Aj iQ B q ] b . 
Thus, assuming the latter condition is met, under the 
action of ideal DD the action of on the code Cj is 



or, in simplified form (assuming /3T m <C 1, ||3>pDD m ||oo < 
1/2, and zero-width pulses): 

F Q [ Ps (T m )Ji^(T m )} (98) 
> 1 - D[p° s (T m ), p^ dcal (T m )] - 2m(cJ 2 + J(3)T 2 . 



proportional to I nj , i.e., is harmless. Quantum logic is 
enacted by the elements of CQ' . Dynamical decoupling- 
operations are enacted via the elements of CQ. We satisfy 
condition (fj[|) because [CQ, CQ'] = 0. 

As an example, consider quantum computation with 
the Heisenberg interaction [H, SH, EH • For the purposes 
of quantum computing with electron spins in quantum 
dots, where a linear system-bath interaction of the form 

H s n B= E E<^> ( 95 ) 

a—x,y,z j 

is the dominant source of decoherence due to hyperfine 
coupling to impurity nuclear spins, it is convenient to 
use only Heisenberg interactions 7?Heis = Si<j Jij ®i ' ®j j 
without physical-level single-qubit gates [42j]. Here &j = 
(<rj, a?, <jj) are the Pauli matrices on the jth system 
qubit and are arbitrary bath operators. To beat 

H-SB we use the Abelian "universal decoupling group" 
Q uni = {I, X, Y, Z}, where X = a*, Y = ®. a), 

Z = a]. It is simple to verify that Ug nni (H^) = 0. 
This is compatible with using Q un i to eliminate Hgg, 
since the global X, Y and Z pulses commute with the 
Heisenberg interaction. I.e., this is an explicit example 
of Eq. ((4]), where we identify -ff c tri with i?Hcis, and i?DD 
with the Hamiltonian generating the global pulses X, Y 
and Z, namely J2j a fi a — x,y,z. As is well known 
[23l |4ll |421 | , universal quantum computation is possible 
using only the Heisenberg interaction provided qubits are 
encoded into appropriate decoherence-free subspaces or 
subsystems. 



VIII. DISCUSSION AND CONCLUSIONS 

Let us now combine our two main results, Eqs. (|56[) 
and (f9"Tj) . for m PDD cycles, each of duration T, i.e., of 
total duration T m , involving N m pulses each of width 5 
and interval r: 

(96) 
(97) 



We remind the reader that the term 
D[p° s (T m ),p^ ea \T m )} is the error due to control 
imperfections in the uncoupled setting, and must be 
dealt with using methods such as fault tolerant quantum 
error correction, composite pulses, or pulse shaping. 



F Q [ Ps (T m ), pT a \T m )] > 1 - D[p%{T m ), P T a \T m )] - min[l, -(e 



211$ 



PDD m || co 



|^pdd,„||oo < c(JT m ) /m + N m J5 + JT m min 

I 



exp(2/3T m /m) - 1 _ 
' 20(T m /m) 
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The term c(JT m ) 2 /m = mcJ 2 T 2 in Eq. {97]) is a bound 
on the error due to the fact that we have terminated the 
Magnus expansion at second order. It can in principle be 
improved by performing a more careful higher order per- 
turbation theory analysis. The term N\ ong JS is the error 
due to finite pulse width. This error can be improved by 
using pulse shaping techniques 0, Hfl ■ The last term in 
Eq. (|97p is a bound on the undecoupled errors, i.e., errors 
due to imperfect decoupling. Considering the zero-width 
pulse limit, Eq. (f9"5)l . we see that provided the number of 

i.e., 



cycles m scales more slowly than [2(cJ 2 + J(3)T 2 
if 

m = o{[2(cJ 2 + Jf3)T 2 }- 1 }, 



(99) 



the fidelity is guaranteed to be dominated by the error 
D[p° s (T m ),p^ eal (T m )} due to control imperfections (the 
"little-o" notation means that the right-hand side domi- 
nates the left-hand side asymptotically). 

We also recall that (3 = ||-ff S cc||oo = ||-ff c tri <8> Ib + 

IS ® H B \\oo < ll-ffctrlHoo + \\H B \\oo and J = ||ffcrr||oc = 

\\Hsb + Hs,res\\oo, where Hsb is the system-bath interac- 
tion Hamiltonian and Hs.rcs are residual undesired pure- 
system terms that do not commute with H ctI \. Express- 
ing the system-bath interaction as Hsb = J2 a <8 B a 
(sum over system times bath operators), we have J < 
J2 a WSaWaoWBaWoa + \\ H s , ICS \\ oc ■ For local Hamiltoni- 
ans involving n system qubits we can reasonably expect 
J oc n (e.g, for electron spin qubits, each of which is 
coupled to a local bath of nuclear spin impurities) . Simi- 
larly, we have ||i? c tri||oo oc n (assuming full parallelism in 
the operation of the quantum computer). The norm of 
the pure-bath Hamiltonian (||i?s||oo) may be very large, 
though in practice it is always finite due to a high-energy 
cutoff or spatial cutoff determining the relevant bath de- 
grees of freedom. Assuming that we are dealing with a 
bath for which IjiJslloo oc Mn (appropriate spatial cutoff, 
such that the n qubits couple to a bath with M degrees 
of freedom, where M can be very large), we also have 
(3 oc n. Thus, we have from Eq. that for fixed T, 



-2+£ 



(100) 



where d is a dimensionless constant involving the various 
energy scales of the problem and e > 0. This last result 
establishes that using PDD with fixed cycle time, there 
is a tradeoff between the number of cycles and the size of 
the quantum register, i.e., there is a limit on scalability. 
On the other hand, the complete inequality suggested by 
Eq. Ijggj is 



mT < [2(cT + J/3)] 



-1/2 



(101) 



so that a better strategy might be to invest resources in 
shrinking the cycle duration T with n, so as to increase 
the number of cycles m. 

Ultimately, based on various comparative studies [ll|, 
HI, HH HH , we expect that there are strategies that will 
outperform PDD altogether and will lead to much im- 
proved scalability. Such strategies are concatenated DD 



0, HI, H EH HI , randomized DD (TSGjllIl, and spe- 
cially tailored DD such as the sequence proposed in [37| 
for the diagonal spin-boson model. We expect that the 
rigorous analysis we have presented here will prove useful 
in the analysis of these more elaborate pulse sequences. 
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APPENDIX 

We provide background and prove the various Lemmas 
found in the main text. For convenience we restate all 
the Lemmas in this appendix. 



APPENDIX A: MAGNUS EXPANSION 

This section is a brief summary of [26|, [43|. The 
Magnus expansion is a method for solving first-order 
operator- valued linear differential equations: 



dU(t,0) 
dt 

17(0) 



= -iH(t)U, <>0, 
= I. 



(Al) 



Here H{t) can be any bounded linear operator. When 
H{t) is hermitian (the only case we consider) Eq. (|A1[) is 
the time-dependent Schrodinger equation and the Mag- 
nus expansion provides a unitary perturbation theory, in 
contrast to the Dyson series. The unitary nature of the 
Magnus expansion is one of it most appealing features. 

The formal solution of Eq. (|A1|) is the time-ordered 
integral 



N 



U(t) 



lim IT 

7^^, J- J- 



exp 



N \N 



T exp[-i [ H(s)ds] 
Jo 



(A2) 



The Magnus expansion represent the solution in the form 
U(t) = exp[— iQ(t)] and expresses ft(t) in a series expan- 
sion. When H(t) commutes with J * H(s)ds the solution 

is U(t) — exp[— i J*H(s)ds}, t > (no time-ordering). 
Otherwise the solution is an infinite series: 



U(t, 0) = lim e 



iM n (t) 



(A3) 
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where M n (t) is the hermitian operator 

n 

M n = *k (A4) 

i=l 

where 

ni[-ff(*)]*o (A5) 

= J2 c j,*[ I ■ ■ I [[H(h),--- ,H(t n )}}dt n ---dh, 
Jti Jt 2 Jt n 



0<*i<-<t„<t 



where [[H(ti),--- ,H(t n )]] denotes an nth level nested 
time-ordered commutator expression between i?(ti), and 
the coefficients Cji are recursively defined and can be 
computed to any order. The first few terms are: 



t 

'o 

ft /■« 



(A6) 



^2 = \ I dt x [ dt 2 [H(t 1 ),H(t 1 )} 
1 Jo Jo 



°3 = / / <K a / dt 3 [H(t 3 ),[H{t 2 ),H(h) 
LZ Jo Jo Jo 



t ft! r t 2 



+ - dtx ds' / dtsimt^Hfa^Hih)]. 
4 Jo Jo Jo 

A sufficient (but not necessary) condition for absolute 
convergence of the Magnus series M n {t) in the interval 
[0,t) is 

f t \\H( s )\\ 00 d s <7T. (A7) 
Jo 



APPENDIX B: EVOLUTION LEMMA FOR A 
SWITCHED HAMILTONIAN 

We prove Lemma [TJ 

The propagator generated by a "switched Hamilto- 
nian" 

H(t) = Hi(t) t^x<t<t u i = l,...,N, 

can be decomposed into corresponding segments: 

U(t N ,t ) = Z7(tjv,tjv--i) • ■ • Z7(*i,to), (Bl) 

where U{b%+\, ti), with tj < t < fj+i, satisfies the 
Schrodinger equation 



dU{t,U) 
dt 



= -iH{t)U{t,U), U(U,ti) = I. (B2) 



Proof. Denote the propagator generated by a time- 
dependent Hamiltonian _ff (t) starting from an initial time 
to by U(t,to). Evolving backward in time from ti to to, 
followed by a forward in time evolution from to to t yields 
a net evolution from t,- to t: 



Letting t = tjv and ti = tjv-i we thus have: 

U(t N , t ) = U(t N , ttf_i)l7(itf_i, to). (B4) 

Repeating this via f7(t/v-i, to) = 
U(tjsr-i,tN-2)U(tN-2,ta) etc. we arrive at Eq. (|B1[) . 
To prove that U(t,U) satisfies Eq. (|B2[) we differentiate 
Eq. (|B3|) with respect to t: 

g/fcji) = dU(t,t ) f 
dt 1 " 0j 
= -iH(t)U(t,t )U(t h t o y 
= -iH(t)U(t,U). (B5) 



APPENDIX C: NORMS AND DISTANCES 

Throughout this work we use unitarily invariant norms 
on bounded operators A [4J|(Ch.4): 



\\A\\ = \\UAV\\ U,V unitary. 



(CI) 



A norm is weakly unitarily invariant if \\A\\ = \\UAU'\\ 
for every unitary U. Obviously, if a norm is unitarily 
invariant then it is also weakly unitarily invariant. In 
addition to being subadditive, i.e., satisfying the triangle 
inequality (by definition of a norm) ||A + P|| < \\A\\ + 

\B\\, unitarily invariant norms are also submultiplicative 

3](p.94): 



Define 



\AB\\<\\A\\\\B\ 



\A\ = VaTa. 



(C2) 



(C3) 



The set of all square matrices, together with a sub- 
multiplicative norm, is an example of a Banach algebra, 
and every C* algebra is a Banach algebra. Note that 
not all matrix norms are submultiplicative. For exam- 
ple, if we define = maxy |a.y| then for the matrices 

J J ) we have ||A|| A = ||£?|| A = 1 but 



A = B = 



\AB\\ A = 2. 

We now give three important examples ;44j(p.91-92). 
The trace norm is 



= tr(|A|) 



trA 



(C4) 



U(t,t t ) = U{t,t )U{ti,t o y . 



(B3) 



Note that if p is a density matrix then ||p||i = trp = 1. 
The trace distance D(pi, p 2 ) = ^\\pi — P2W1, plays a spe- 
cial role since it captures the measurable distance be- 
tween different density matrices pi and p 2 [4a] . Namely, 
D(pi,p 2 ) is an achievable upper bound on the trace dis- 
tance between probability distributions arising from mea- 
surements P performed on p\ and p 2 [l[ (Theorem 9.1), 
in the sense that D(pi,p 2 ) = ma.x P ((P)i — (P) 2 ), where 
P < I is a positive operator, and (P)i = tr(PpA 
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The Frobenius (or Hilbert- Schmidt) norm 



First note that 



\\A\\ 2 = x /p^4) = v /tr^4U) 

= (Em 2 ) i/2 > 



(C5) 



(where A has matrix elements a^) is the norm induced 
by the Hilbert- Schmidt inner product 



{A, B) = tvA^B. 



(C6) 



Finally, the operator norm is ||^4||oo = s i(A) = 
su P||^||=i II^IV*)!!) where si(A) is the first (largest) sin- 
gular value of A, i.e., the largest eigenvalue of \A\. 

All three norms are multiplicative with respect to the 
tensor product Hi|(Ch.2): 

||i4<8fl||i = HAHillBlli i = l,2,oo. (C7) 

They satisfy the ordering 

imioo < H^iia < imii. (cs) 

Another useful inequality is [46((Ch.2) 

\\ABC\\ < IIAIIoollBlHICIIoc, (C9) 

where ||-|| denotes any unitarily invariant norm. A special 
case of this is obtained by setting A, B or C = I: 

\\AB\\ < WAIUIBIWAW^WBW^WBIUIAW. (CIO) 

An important inequality we need relates the norm of the 
partial trace and the norm of the operator being traced 
over (for a proof see [3ll|): 

||tr B A|| 4 < di\\X\\i (i = l,2,oo), 
di = 1, 



d 2 = ^dim(7i B ), 
d x = dim(7Y B ), 



(Cll) 



where A is a linear operator over the tensor product 
Hilbert space Hs^Ti-B- For the trace norm this is a spe- 
cial case of the well known result that trace-preserving 
maps (in this case the partial trace) are contractive [l[. 



APPENDIX D: NORM TO ERROR PHASE 
INEQUALITY FOR MIXED STATES 

We prove Lemma [21 

Let U = exp(— iA) where A is hermitian. Then for any 
sub-multiplicative norm 

\\UBU f -B\\ < \\B\\ min[2,e 2||A||o ° - 1] 

2||A|U<1 

< 2||B||min[l,(e-l)||A|| 00 ] 

Proof. 



e x - 1 = x 1 + 



X<1 

< X 



(e-l)x. (Dl) 



C<1 



By a similar calculation we also get e — 1 < (e — 2)x. 
By the triangle inequality 

\\UBU* - B)\\ < \\UBU*\\ + \\B\\ < 2\\B\\. (D2) 

On the other hand, using the Taylor expansion of 
cxp(— i[A, •]) we have 



\UBW - B\ 



n=0 



[ n A, B] - B\ 



= iiE^^n 

< Y.h\[ n A,B]\\ 



< 



E 



\B\\ 



n=l 

( e 2||A|U _ jjpn 

2||^||=o<l 

< 2(e- 1)||^|| oo||S||, (D3) 

where in the penultimate inequality we iterated 

\\[A,B}\\<2\\AB\\<2\\A\\ 00 \\B\\, (D4) 

[where we used submultiplicativity together with 
Eq. fCTO]) ] to get 



[ n A > B}\\<\\A\\ 00 \\[ n ^A,B}\\<...<\\A\U\B\ 



(D5) 



APPENDIX E: MAGNUS EXPANSION 
TRUNCATION BOUND 

We prove Lemma [U 

Consider a time-dependent Hamiltonian H(t), < t < 
T, and the partial sum of fcth and higher order terms in 
the corresponding Magnus expansion: 



*fc = E ^ 



Assume the Magnus expansion converges in the trace 
norm. Then 

||$*|| <c k (T sup ||i?W||oo) fe 

0<t<T 
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where Ck = 0(1) is a constant. 

Proof. Define h = sup 0<t<T H-H^t)!^ and rescale H{t) 
by hT: 



H\t) = 



H(t) 
hT ' 



(El) 



We can rewrite J7, as: 



Jo Jo 



[[H'(h)> ■ ■ ■ , H'(t n )]]dt n ■ ■ ■ dt 1 
= (hT) l n % lH'(t)]T. (E2) 

Recall the condition for absolute convergence of the 
Magnus expansion, Eq. (fAT)) . Since J Q T \\H(t)\\ dt < 
T sup 0<t< y ||J?(t)|| a sufficient condition is: 



hT < 1. 



(E3) 



Absolute convergence (convergence of the sum of abso- 
lute values) means that if we define, for k > 1, the partial 
sum 



B„, fc ^ || ft, [#(*)]£ 



(E4) 



then 

lim B n>k =p k [H{t)]l < oo, (E5) 

n— *oa 

where (3k is some functional of H(t). Similarly for H f (t): 

n 

lim Y, lW(*)]o II = 0k[H'(t)]o = ^ = 0{l). 



i—k 



/E6) 

Let us now focus on the partial sum of fcth and higher 
order terms in the Magnus expansion, $fc = ^2 i=k ^i- 
We can bound in the following manner: 



i* fc ii < E w^mm 



■i—k 



Eq.JH 



(hT) k Y(hTy- k \MH'(t))T 



Eq.JH «E1 

< (fcr)*x;ii n *[^ , (*)]oi 

i=fc 



(E7) 
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